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Abstract 

Estimation of extreme value copulas is often required in situations 
where available data are sparse. Parametric methods may then be the 
preferred approach. A possible way of defining parametric families 
that are simple and, at the same time, cover a large variety of mul- 
tivariate extremal dependence structures is to build models based on 
spectral measures. This approach is considered here. Parametric fam- 
ilies of spectral measures are defined as convex hulls of suitable basis 
elements, and parameters are estimated by projecting an initial non- 
parametric estimator on these finite-dimensional spaces. Asymptotic 
. ^ distributions are derived for the estimated parameters and the result- 

^ ing estimates of the spectral measure and the extreme value copula. 

Finite sample properties are illustrated by a simulation study. 



1 Introduction 

Extreme value copulas provide a suitable general approach to modelling mul- 
tivariate extremes. Various nonparametric methods for estimating extreme 
value copulas have been proposed in the last few years [T7], [10], [H], [3] 
(also see [15], [5] and [13] for related approaches). In practical applications, 
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such as for instance operational risk or rare natural disasters, one is however 
often in a situation where available data are sparse. Nonparametric methods 
generally require a fairly large sample size in order to be reliable. For small 
samples and in situations where one may have some idea about plausible 
properties of the distribution, parametric methods are likely to yield more 
accurate results. An approach to parametric inference for extreme value 
copulas is discussed for instance in [2]. 

One of the key issues is how to define parametric families that are simple 
and at the same time general enough to cover a large variety of multivariate 
dependence structures in the extremes. For instance, some of the most pop- 
ular models are based on Archimedean copulas, which all correspond to the 
same type of extremal dependence structure, characterized by the Gumbel 
copula [JJ. One way of achieving more flexibility in the extremes is to build 
models based on spectral measures. This is the approach taken here. For 
related work see e.g. [6], [12], and [TT] . 

More specifically, the idea pursued in the following is to select a finite 
number of suitable spectral measures as basis elements and to use their con- 
vex combinations as a parametric family of dependence structures. Given a 
sufficiently large number of such basis elements, any spectral measure can be 
approximated by a weighted sum. Estimation of the coefficients can then be 
carried out by projecting a nonparametric estimator, such as the one in [3J, on 
the finite-dimensional space generated by the basis elements. If the number 
of basis elements in the model is large (and increasing with the sample size), 
then projecting the original non-parametric estimator can be considered as 
a discretization technique. This is the setting in (6] and 

On the other hand, an appropriate model with a small number of basis 
elements can have the advantage of dimension reduction. Given a reason- 
able parametric model with a small number of parameters, one can reduce the 
variability of a nonparametric estimator by projecting it on a low-dimensional 
space. This is the approach studied here. We define explicit parameter esti- 
mators in the low-dimensional setting and study the asymptotic distribution 
of the resulting estimators of the dependence structure. To illustrate the 
potential advantage of dimension reduction, we construct an example with 
three basis elements and compare a non-parametric estimator with its low- 
dimensional projection in a simulation study. 

Note that in principle any nonparametric estimator (cf. [5], [3], [13], [TT] . 
|lUj . [TT] ) can be used as a starting point. Depending on the nonparametric 
method used in the projection, the marginal distributions are either known 
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or estimated from the observed data. The asymptotic results given below 
only require that a functional limit theorem in a suitable topology holds for 
the initial estimator. 

The paper is organized as follows. Basic definitions and concepts of multi- 
variate extreme value theory are summarized in section 2. Parametric models 
in the spectral domain and a corresponding parametric estimator are intro- 
duced in section 3. Asymptotic results, including consistency and a central 
limit theorem, are derived in section 4. The theoretical results are illustrated 
by simulations for a specific model in section 5. Final remarks in section 6 
with a discussion of some open problems conclude the paper. 

2 Basic definitions 

Consider a sample X 1; . . . , X n consisting of iid realizations Xj = (Xj i, . . . , X itd ) 
of a rf-dimensional random vector X = (Xi, . . . ,X d ) T G M. d with marginal 
distributions F 1 ,...,F d and copula Cx- That is, Fj(t) = P (Xj <t) for 
j — 1, . . . , d and t G R, and 

P(X<x) = Cx(F 1 (x 1 ),...,F d (x d )) 

for x G M. d . The notation x < y for x, y G R. d means Xj < \jj for j — 1, . . . , d. 
The transposition operator (-) T in X = (Xi, . . . ,X d ) T indicates that X is 
considered as a column vector. Distinguishing columns and rows will be 
useful in some calculations later on. 

The vector M n = (M n l , . . . , M n d ) T of componentwise maxima 

Mr, j = max Xj ~ 

n ' J i=l,2,...,n ,J 

then has marginal distributions P(M n j < t) — F- l (t) and a copula C Mn (u) 
given by 

P (M n < x) = C Mn (F? ( Xl ) ,...,F% (x d )) = CI (F, ( Xl ) , . . . , F d (x d )) . 

In the limit one obtains, under general conditions, an extreme value copula 
C (u) = lim^oo C x (u 1 /™) (u = («!,..., u d ) T G [0, l] d ) with the characteris- 
tic max-stable property 

C (u) = C n (u 1/n ) (1) 
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for all n G N. For an accessible introduction to this topic see e.g. [HJ and 
references therein. The definition ([I]) of extreme value or max-stable copulas 
is equivalent to the representation 

C (u) = exp (-^(-logwi, • • ■ , -logttd)) 
with the tail dependence function 

£(x) = / max (tu^) {w 1 , . . . , w d ) , (x G [0, oo) d ), (2) 
J Ad i=i,...,d 

and \1/ the so-called spectral measure on the unit simplex in M d , A^ = {x G 
[°> !] d : E?=i ^ = !}> satisfying 

/ Wid^f (wx, . . . ,Wd) = 1, (i = l, (3) 

^A d 

(cf. Theorem 6.2.2 in [9]). Note that the last condition implies f A d^fr (w) = 
d. The underlying original results go back to [I] and [To]. The main conclu- 
sion is that spectral measures, tail dependence functions, and extreme value 
copulas are equivalent representations of dependence structures in multivari- 
ate extreme value theory. Note also that £(rx) = r£(x) for r > and 
x G [0, oo) d , so that it is sufficient to specify £(x) for x G A^ only. The re- 
striction of i to Arf is also called Pickands dependence function and is usually 
denoted by A (•). The extension to x G [0, oo) d is obtained by 

*(x) = Uxll^pj-) = \\*\\iA(w) (4) 

where ||x||i = x\ + ■ ■ • + x<i (since Xj > 0), and w = x/||x||i. Note that 
condition (|3| is equivalent to 

A(ej) = l, (j = l,...,d) (5) 

where = (e^i, ...,ej^) T is the j-th unit vector in R d : ejj = (j ^ I) and 
Cjj = 1. That is, ^ standardizes A on the vertices of the unit simplex A^. 
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3 Parametric models for spectral measures: 
construction and estimation 




3.1 Models 

One way of building parametric models that encompass a large variety of 
extremal dependence structures is to start at the level of the spectral measure 
Thus, suppose that ^f%. . . ,^ p are some fixed spectral measures. The 
correspoding dependence functions and extreme value copulas will be denoted 
by £\, . . . ,£ p , Ai, . . . , A p and C%, . . . , C p , respectively A parametric family of 
spectral measures V p = (-,0) ,9 E Q}, and corresponding families A p and 
C p of (Pickands) dependence functions and copulas respectively, can then be 
obtained by defining spectral measures of the form 

p-i / 

*o,0) = 5>*<(o+ i 

where 9 = (0 X , . . . , G 6 and 6 = {tf G (0, l)*" 1 : Ef=i ^ < !}■ As 
^ remains valid for convex combinations, \I/ (•, 9) is a spectral measure by 
definition. In terms of the corresponding dependence functions we have 

p-i 

e(x,9) = J2 6 ^ W- 

1=1 

p-1 

A(w,0) = £Mi (w) 

For the copulas we obtain 

p-i 

C (u, 0) = exp <( - V Ot % (- log u) - l p (- log u)] - (- log u) J. (6) 





where ^ (-logu) = t { (-log«i, . . . , - log« d ). 

Henceforth we assume that the parameter # is identifiable in the sense 
that Ag = Agi implies 9 = 9' for 9,9' G 0. A sufficient criterion for the 
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identifiability of 9 is linear independence of the basis elements Ai, . . . , A p . If 
9 is not identifiable, then an estimator 9 may fail to converge. An important 
example of this issue is the decomposition of discrete spectral measures. 
According to [H] , any discrete spectral measure on A2 can be expressed as a 
convex combination of two-point spectral measures. This result can also be 
written in terms of piecewise linear dependence functions and Marshall-Olkin 
copulas. However, the decomposition is not necessarily unique. This can be 
illustrated by the following example. Let (t, 1 — t) G A 2 be represented by 
the first coordinate t G [0, 1] and consider the family Q 4 of discrete spectral 
measures \1/ = Y^t=i c «^(i-i)/3 with c$ > for % = 1, . . . ,4. It is easy to see 
that a basis of 2-point spectral measures needed for the decomposition of all 
f G Q 4 must include all elements of Q4 with only two atoms. These are 

13 3 1 

^1 = S + Si, #2 = 6y 3 + 82/3, ^3 = + 2^2/3, ^4 = 2^1/3 + 2^1- 

The non-uniqueness follows from + = ^(^3 + ^4). 
3.2 Estimation 

Several nonparametric estimators of Pickands dependence functions, spec- 
tral measures, and corresponding extreme value copulas have been proposed 
in the recent literature [5], [3], [13], [17], [10], [11]. Generally, these meth- 
ods require fairly large sample sizes in order to achieve a sufficient degree of 
accuracy. In contrast, parametric estimates are expected to be reasonably 
accurate for moderate or even small sample sizes, provided that the paramet- 
ric assumptions are sufficiently realistic. To see how much may be gained by 
parametric estimation, we consider the following approach. Suppose that a 
nonparametric estimate A of the dependence function A is given, and recall 
that the corresponding I is obtained from A according to Q. A natural 
parametric estimator based on the family A p , and A as intial estimate, is 
obtained by projecting the function A on A p . Note that even if @ does not 
hold for A (see e.g. [10J ) , it holds automatically for the projection of A on A p , 
so that this projection is a proper dependence function by definition. Note 
also that the projection improves the accuracy of the estimate if the true 
dependence function is indeed in A p . Related improvements for projections 
on infinite-dimensional spaces of spectral measures and approximations by 
sieve methods have been considered in [6]. 
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Specifically, we may start for instance with the following nonparametric 
estimator A considered in [10], p2], [3], [5] and [15]. Suppose that the depen- 
dence structure of X G M. d (i = 1, . . . , n) is characterized by an extreme value 
copula C (-,9) G C p , and, as before, the marginals are denoted by Fx, ... , F d . 
Given n iid realizations Xj (i = 1, 2, . . . , n), define Yj = (Y it i, . . . , Y i)d ) by 

Ytj = - log FjiXij). 

Then Yy are standard exponential random variables, and 

P (Y hl > Vl ,..., Y i4 >y d ) = C (e"* e~^,9) = exv(-£(y, 9)) 

= exp(-||y|| 1 A(w,0)) 

for y G [0, oo) d (and w = y/ 1 1 y 1 1 x G A^). Hence, for 

Yi ■ 

£«(w) := min — — (i — l,...,n) 

j=l,..,d Wj 

with w G Arf one obtains 

P(&(w) > t) = P (y ijX > Wl t, . . . , Y itd > w d t) = exp(-M(w)). 

This means that, for any fixed w G A^, £i(w), . . . , £ n (w) are iid exponen- 
tially distributed with mean 1/A(w), and — log(^(w)) (i = 1, . . . , n) are iid 
Gumbel distributed with location parameter logA(w). In particular, 

E(-log6(w)) =logA(w)+ 7 

where 7 is the expectation of the standard Gumbel distribution (i.e. 7 = 
r'(l) ~ 0.5772, the Euler-Mascheroni constant). A nonparametric estimator 
of A may therefore be defined by (see [15], [5], [T3], [6], [TO] ) 

1 n 

logA(w) = — ^log&H-7(wG A d ). (7) 

i=i 

Unfortunately, A satisfies §2§ and ([5]) only by chance. The standardiza- 
tion ([5]) can be achieved by modifications of log A that substract a suit- 
able linear combination of logA(w) evaluated at certain values of w (cf. [3], 
[17] . [TU]). In particular, [TU] defines a nonparametric least squares estima- 
tor (nonparametric OLS) Aols(w) by log Aols(w) = /3 (w) where /3q(w) 
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is obtained by least squares regression of log^(w) — 7 (z = 1,2, ...,n) on 
— log£j(ei) — 7, — logged) — 7. The resulting estimate satisfies (|5]), but 
it still may fail to satisfy ([2]). 

In contrast, in the parametric approach introduced above a modification 
is not needed, because the projection on A p automatically leads to a proper 
Pickands dependence function. However, the examples below demonstrate 
that the initial estimator remains crucial for both the asymptotic distribution 
and the finite sample behaviour of the parametric projection. Specifically, 
we will compare the parametric approach based on A and Aols respectively. 

Classical results on empirical processes yield a functional central limit 
theorem of the following form (see [TU] and references therein). Let C(A^) 
denote the Banach space of real- valued continuous functions on A^ equipped 
with the supremum norm. Then, as n — > 00, 

V™ (A (w) - A (w) j A A (w) ( (w) =: CA,non P (w) in C(A d ) (8) 

where ( (w) (w G A^) is a zero mean Gaussian process with covariance 
function 

7? (v,w) := cov(C(v),C(w)) = cov(-logf(v),-logf(w)) (v,w G A d ). 

Note that the joint distribution of £i(v),£j(w) does not depend on i, so that 
dropping the index i here does not lead to confusion. This result implies 
in particular that, for large n and Wi, . . . , wjy G A^, the joint distribution 
of A(wi), . . . , A(wn) can be approximated by an iV— dimensional normal 
distribution with mean 

fl (w 1; . . . , Wat) = (A(wi), . . . , A(w N )) T 

and covariance matrix = n~ l \o (\Vj, w J -)] i =1 N with 

a (Wi,Wj) = A(w i )A(w i )7c ( w ^ w j) 

= ^(wO^WjOcovC-bg^wO.-log^WjO). (9) 

Now, given a parametric class of spectral measures V p based on ^1 ... , ^/ p , an 
estimator of the corresponding Pickands dependence function A(-,9)e A p 
can be defined as follows. Let A be the preliminary estimator in (M), and 
denote by Ai (i — 1, . . . ,p) the Pickands dependence functions correspond- 
ing to the spectral measures {i = 1, . . . ,p). Define a grid of w- values 
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w°, ...,w% £ A d , and the vector a = (a®, . . . , a° N ) T with a° = a(w°) 
(i = 1, . . . , N) and 

a(w) = 1 (w) - A p (w) , w G A d . 

Furthermore, define the iV x (p — 1) matrix H° = [hj(w^)] i=1 N . j=1 p _ 1 
with 

h j (w) = A j (w)-A p (w), w G A d . 
Then the least squares estimator of 6 1 is equal to 

9 = g°a° 

where Q° = (H° T H )' 1 H° T . Since the dependence functions Aj (w) are 
defined for all w G A d , an estimate of A (w) is available for any w G A p by 
setting 

A(w, 0) = A p (w) + h T (w) 9 = A p (w) + h T (w) Q°a° 

where h (w) = (hi(w), . . . , h p _ 1 (w)) T . 

Letting N tend to infinity, an estimator of 9 based on all values in A d 
can be obtained as follows. Suppose that the grid wj,...,w^ is chosen 
such that, as N — > oo, the point measure Mn(B) = N" 1 J2f =1 1 {wf G B} 
(B G B(A d )) converges weakly to a probability measure M on A d with 
Lebesgue density m(-). This can be achieved by deterministic or by random 
choice (by sampling from M) of the grid. Then 

N . 

N-Uh^H ) =N- 1 J2Hrf)h J (rf)= /ii(w)^(w)dMjv(w) 



converges to 



s i,j 



= / hi(w)hj(w)m(w)dw. 



This follows from the continuity of all Aj (and hence all hj) and from the 
compactness of A d . Similarly, the limit of 

N 

3 i=i 
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IS 



r .i 



hj(w)a(w)m(w)dw. 
Thus we obtain an estimator that depends on the density function m, 

Om = S-\ 

where S = (si 1 j) i - 1 p _ l and r = (n, . . . , r p _i) T . Even more generally, the 
previous estimators can be seen as special cases of 

Om = S-h (10) 

where 

Si.j = / hi(w)hj(w)dM(w), 
hj(w)a(w)dM(w) 



and M is any distribution function on such that S is of full rank. 

The same approach can be applied to any initial nonparametric estimator 
of A for which a functional limit theorem is available. In particular, for the 
nonparametric OLS, A OLS , Gudendorf and Segers ([10]) obtain 

(Aols (w) - A (w)) A A (w) [( (w) - A^(w)C (e)] 

= : CA,OL5,nonp(w) in C(A d ) 

where ( (w) is the Gaussian process defined in m), £ (e) = (£ (ei) ,...,( (ed)) T , 
A opt (w) = Yr x E [C (e) C(w)] and £ = E [( (e) C (e)] . Applying the paramet- 
ric approach, we define as before 

9m,ols — S 1 ?ols (11) 

with 



and 

a<3Ls(w) = Aqls (w) - A p (w) , w G A d . 
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4 Asymptotic results 



We will use the notation er (v, w) = cov (C4,non P (v), Cvonp(w)) and a LS (v, w) 
cov (CA,OLS,nonp( v )) CA,OLS , ,non P ( w )) for the asymptotic covariance functions of 
A and A OLS respectively. The asymptotic distributions of 9 m and 9m,ols are 
given by 

Theorem 1. Let X; = (X^, . . . , X i4 ) T e R d (i = 1,2, ... ,n) be iid re- 
alizations of a d- dimensional random vector X with marginal distributions 
Fi,...,Fd and extreme value copula C(-,9°) G C p . Denote by Aj (j = 
1, . . . ,p) the Pickands dependence functions defining C p , and let 9 m be defined 
by ilU) and 9m,ols by (11), where M is such that S is of full rank. Suppose 



furthermore that the parameter 9 is identifiable and 9° is in the interior of 
the parameter space O = j# = (9±, . . . ,9 p _i) T 6 : H^j^ < 1 \. Then, as 

n — > oo, 9m and 9m,ols converge to 9° in probability, and 



4Z, 



^ [8m,ols — #°J - > Zols 

where Z and Tiqls are {p ~ I) -dimensional normal random vectors with zero 
mean and covariance matrices 

cov(Z) = V = S^tt {S^f , 
cov(Z 0i5 ) = Vols = S^Qols {S^f 
where Q = [^j,i]j,i=i,..., p -i and ^ols = j i=i ...,p-i are de fi ned b U 

u>j,r.= [ [ h j {v)h l {w)a(v,w)dM{v)dM{w), (12) 

hj(v)hi{w)a L8(-v, w) dM{v)dM{w). 



OLS ._ 



Proof. Since the proof for 9m and 9m,ols is the same, it is stated for the first 
estimator only. We have 

9 M = S' 1 J h(w) - A p (w)) dM(w) 

= S- 1 <f> (A-A^J, 
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where 4>(f) = f A h(w)/(w)dM(w) is a linear mapping from C(A^) into 
R p_1 . Analogously, we have 8° = S ,_1 0(v4 — A p ), and hence 



y^^-^oj = S- X (j) (A -A^. 



It is obvious that the mapping / i-> S~ 1 (p(f) is continuous. Hence the func- 
tional Central Limit Theorem ([8]) for A and the Continuous Mapping Theo- 
rem yield 

Vn(0-0°) A,S-V(Ca ot ) =: Z. 

Recall that C4,nonp is a zero-mean Gaussian process with covariance func- 
tion cr(v, w) = E[(A,nonp(v)(A,nonp{"w)) introduced in Hence, as a linear 
mapping of (A.nonp, the random vector Z is Gaussian with zero mean and 
covariance matrix V = 5 ,_1 r2(S' _1 ) T , where Q = [Wy]tj=i,...,i-p is the co- 
variance matrix of <p((A,nonp) ■ The representation (12) follows from Fubini's 
Theorem: 

= £ / /ii(v)CA,tump(v)tiAf(v) / /i i (w)C j4inonp (w)rfM(w) 

jA d JA d 
= 11 ^(v)/l i (w)E[CA,™ m p(v)CA,no n p(w)]dM(v)dM(w). 

□ 

An immediate consequence of this result is the asymptotic normality of 
A(w, 9 M ) = A p (w) + h T (w) 9 M and A(w, 9 m ,ols) = A p (w) + h T (w) 6 m ,ols 
uniformly in w G A^. 

Corollary 2. Under the assumptions of Theorem^ we have, as n — » oo, 

V^(a(w,9m) -A(w,fl )) 4(,(w), (13) 

(A (w, §m,ols) - ^ (w, 6"°) J A Ca,ols (w) 

where (a (w), Ca.ols ( w ) T w are zero-mean Gaussian process with co- 

variance functions 

7a (v, w) = h T (v) Vh (w) , 

and 

7a (v, w) = h T (v) VoLsh (w) , 
where V and Vols ore as in Theorem [7} 
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Proof. Recall that 



A(w, 9 M ) = A p (w) + h 1 (w 



'M- 



As the mapping 9 H- h T 9 is linear and continuous in C(Ad), we obtain (13) 
from the Continuous Mapping Theorem. In fact, we have the representation 
O(w) — h T (w)Z in (13). The covariance structure of the limit process 
follows from 

cov (h T (v)§ M ,h T (w)9 M ^ = h T (v)var(^ M )h(w). 

□ 

Note that, more specifically, Theorem [l] implies that y/n(A(w, §m) — 
A(w,9 )) and y/n(A(w, 9m,ols) — ^4( w ;# )) are asymptotically equivalent 
to the stochastic processes Ca( w ) = h T (w)Z and (a,ols( w ) — h T (w)Z Ls, 
respectively, with index w 6 A<j. The random variables Z and Zols are the 
weak limits in Theorem Q] 

Another consequence of Theorem [T] is the asymptotic distribution of 
C(u,9 M ) and C(u,9 M:OLS ). 

Corollary 3. Under the assumptions of Theorem^ we have, as n — >■ oo, 
V^(c(u,9 M )-C(u,9 )) ACc, 

(C(u, 9 m ,ols) - C(u, fl )) A Cc ,OLS 
where (c, (c,ols ore zero mean Gaussian processes with covariance functions 

7c(u,v) = C>,0 o )V(7M ), 
7c,OL5(u,v) = C t (u,9°)VolsC(v,9°). 

Here, 

and V, Vols are as in Theorem^ More specifically, 

Cc(-) = C T (;0°)Z, (c,ols(-) = C T (;9°)Z OLS 
with Z, Zqls from Theorem^ 
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Proof. Recall ^ and denote r = r(u) := || — logu||i and w = w(u) := 
r _1 (— logu), where log(u) is understood componentwise. Then we obtain 

C(u, 9) = exp (-r (fl T h(w) + A p (w))) , 

and hence, for j, k — 1, . . . ,p — 1, 

dg.C(u,e) = -rh ] (w)exp(-r(9 T h(w) + A p (w))) 
d ej e k C(u,6) = r 2 h j (w)h k (w)ex V (-r(9 T h(w)+A p (w))). 

A well known consequence of ^ is that each Pickands dependence func- 
tion A assumes values in [1/d, 1] only. Since 9 T h(w) + ^4 p (w) = A(w, 9) is a 
proper Pickands dependence function, and each hj is a difference of two, we 
obtain 

\d Bj C(u,9)\ < 2r(u) exp(-r(u)/d) 
\d 6] e k C(u,9)\ < 4r 2 (u)exp(-r(u)/d). 

Thus the first- and second-order derivatives of C(u, 9) with respect to 9 
are uniformly bounded in u G [0, l] d , and the Taylor approximation 



C(u, (9) - C(u, 6°) = C(u, 9°)(9 j - 9]) (c (u, 9°^ (6 -9 ) + O( 



9-9 



2 





with C(u, 9) := (d 6l C(u, 9),..., d ep _ x C{vL, 9)) T is uniform in u e [0, 1]. The 
final result now easily follows from Theorem [Tj □ 



5 Examples and simulations 

5.1 A parametric model example 

To illustrate how one may construct spectral measures . . . , ty p for building 
parametric models, we consider an example with d = 2 and p = 3. For ease 
of notation, we parametrize A2 by the first coordinate, so that t G [0, 1] 
represents (t, 1 — t) G A 2 , and we can write d^/(t), A(t), etc. In particular, if 
a spectral measure \1/ has a Lebesgue density /, then the standardization (j3]) 
reads as 

1= / t/(t)dfc= / (l-t)f(t)dt. (14) 
Jo Jo 
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Let 



fi(t) 



|(l-cos(37rt)) xe [0,2/3] 



ab(l + cos(3vrt/2)) t e (2/3, 1] 

/a(t) = /i(l-t) 
f 3 (t) = csin(7rt). 



With appropriate constants a, b, and c, the functions /$ satisfy (14), and 



we define the basis elements of the parametric family by (Nfi(t) = fi(t)dt. 
Figure [I] shows plots of the spectral densities /», fc, fs- The constants a, b, c 
are derived as follows. 



Spectral densities f; 




Figure 1: Spectral densities fi, f2, and fe. 



Due to symmetry, we have Ltfz{t)dt = J^(l — t)fs(t)dt 7 so that (14) 
yields 

f 1 f 1 2c 

2 = / Mt)dt = c sin (irt) dt = —. 
Jo Jo 7T 

Thus, C = 7T. 

To determine a and b in fi and fi-, it suffices to consider j\. Let <7i(z) := 
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f*tfi(t)dt. For z G [0,2/3] one has 



't 2 ' 


z 


t sin(37rt) 


z 


cos[37rt] 


1 


_2_ 





3tt 





9tt 2 


0/ 



a ( z 2 2;sin(37T2;) cos(37rz) 1 
2 V ~2~ 3tt 9tt 2 + 9^2 



and for z e (2/3,1], 



01 (*0 = g + ab 



+ a6 



+ 



J 2/3 



tsin(37rt/2) 
3772 



+ 



J 2/3 



cos(3tt£/2) 
9tt 2 /4 



2/3 . 



2 2; sin(37rz/2) cos(37rz/2) 



9 



+ 



3tt/2 



+ 



1 



In particular, 



a , / 5 2 

= - + abi h , , . 

y K J 9 V 18 3tt 9tt 2 7 



9tt 2 /4 9tt 2 /4 
4 \ 



Moreover, note that §*{l — i)f\{t)dt = h\{z) —gi(z) with hi(z) := J Z f\{t)dt. 
For z6 [0,2/3] one obtains 



hAz) 



z — 



sin(37rz) 



3tt 



and for z e (2/3,1], 

fci(z) 



a6 



t + 



sin(37rt/2) 



3vr/2 



J 2/3 



abi z -2/3 + 



In particular, 



^(1) 



sin(37rz/2) 
3772 

a(l + 6(1 - 2/tt)) 



We need a and 6 such that fi satisfies (14), which is equivalent to <7i(l) = 1 
and h\(l) = 2. The latter equation yields a = 6 (1 + 6(1 — 2/ir))~ 1 . Substi- 
tuting this in <7i(l) = 1, 



7T 



8 - 6tt + 2tt 2 
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and hence 

12?r 2 - 36vr + 4£ 



3tt 2 - 8tt + 8 

Note that further spectral densities of this type can be defined, for instance, 
by replacing 2/3 in the definition of f\ by other values (in the interval (0, 1)). 



5.2 Sampling technique 

The asymptotic results obtained in section 4 are illustrated by simulations 
for the example introduced above. Thus, the copula C(-, 9) is defined by the 
Pickands dependence function 

A(t, 9) = 9 1 A 1 (t) + 9 2 A 2 (t) + (1-9,- 9 2 )A 3 (t) 

with t G [0, 1] representing (t, 1 - 1) G A 2 , < 9t, 9 2 < 1, Bi + 9 2 < 1, and A; 
i — 1, 2, 3 being the Pickands dependence functions corresponding to ^ (and 
/i). A random vector (Xx,X%) ~ C(-,9) can be simulated exactly using the 
algorithm proposed in [5|. More specifically, given a bivariate dependence 
function A = A(-,9), the corresponding extreme value copula C = C(-,9) 
can be sampled as follows: 



1. Simulate Z G [0, 1] with distribution function 

P(Z<z) = z + z(l-z)^=: G z {z). 

Note that if Z has a density gz, then gz = G' z . 

2. Calculate 

Z(l - Z)A"(Z) 
Pl j A(Z)G^(Z) • 

Let V = Ui with probability p(Z) and = U1II2 with probability 
1 — p(Z), where Ui, U 2 are independent and uniformly distributed on 
[0,1]. 

3. SetXi = 1/ Z / A ( Z ) andX 2 = Then the distribution function 
of the random vector (Xi,X 2 ) is equal to C. 
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The computation of Gz, G' z , and p(Z) can be simplified as follows. Recall 
that, given a spectral density / = + #2/2 + (1 — Oi — 62) fz, ^.(2) = 
/ max(fz, (1 - t)(l - z))f(t)dt. Since tz > (1 - i)(l — z) is equivalent to 
£ > 1 — z, we obtain 

A(z) = / (1 - z)(l -t)f(t)dt+ / ztf(t)dt 

JO Jl-z 

= (1 - z)(/i(l - z) - - z)) + z(l - g(l - z)) 

= Z-g(l- Z ) + (l-z)h(l-z), (15) 



where g(z) = J* tf(t)dt and h(z) = f* f(t)dt (note that (14) implies (7(1) = 1 
and h(l) = 2). From (15) we obtain that 

A'(z) = 1 + g'(l -z)-{l- z)h'{l -z)- h(l - z). 

Consequently, g\z) = zf(z) and h'(z) = f(z) imply that 

A'(z) = 1 + (1 - z)f(l - z) - (1 - z)f(l -z)- h(l - z) 

= l-h(l-z), 
A"(z) = h'(l-z) = f(l-z). 

Thus we obtain 

_ zA(z) - (1 - z)zA'(z) _ z(l - g(l - zj) 
z{ ) ~ W) " W) ' 

and therefore 

l- g (l-z) + z(l-z)f(l-z) zjl - g{\ - z)){\ - h{\ - z)) 
Gz[Z) ~ A(z) " A\z) 

It is obvious that all functions /, A, g, h corresponding to C(-, 9) are con- 
vex combinations of the corresponding f i} A^g^hi with weights 0i,9 2 , and 
1 — 9\ — 9 2 . Since /(z) (z G [0, 1]) is a convex combination of the bounded 
functions fa (i = 1,2,3), we can simulate Gz by rejection sampling. Con- 
sequently, for the simulation of (Yi,Y 2 ) ~ G(-,9) we only need 9 and the 
functions fi, A iy g^ hi for i — 1, 2, 3. The representations of /1, g\, and foi are 
already derived above, in subsection |5.1| An explicit representation for A\ 
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follows from (15). Due to f 2 {z) = f\{l — z) one obtains A 2 (z) = A 1 (l — z). 
Furthermore, 

h 2 {z) = [ h(l-t)dt= f f 1 {y)dy = h 1 (l)-h 1 {l-z) 

JO Jl-z 

= 2-h l {l-z), 

92(z) = I (1 - y)fx(y)dy = h^l) - h(l - z) - ( 9l (l) - 9l (l - z)) 

Jl-z 

= l-h l {l-z)+g 1 {l-z). 
Finally, for i = 3 we have 

hs(z) = / fs(x)dx = 1 — COs(7T2) 





CZ 

9s( z ) = I x fs( x )dx = — sm(rrz) — zcos(ttz). 



7T 



Figures |2pi) through d) show typical samples Xj = (Xii,X i2 ) ~ C(-,6) 
(i = 1,2, ...,n) with n = 1000, and 6 = (1,0), (0, 1), (0,0) and (|, |) respec- 
tively. Image (and contour) plots of two-dimensional kernel density estimates 
for these samples are shown in figures^) through 2d). 

5.3 Simulation results 

To study the finite sample performance of the estimators of A discussed 
above, the following simulation study was carried out. For 9 = (0.1,0.1), 
(0.05,0.9) and (0.8,0.1) respectively, 1000 simulated samples of size n = 
25 • 2 J (j = 0,1,..., 8) were generated. For each sample, the nonparamet- 
ric estimates A and Aols as well as the corresponding parametric esti- 
mates A(w,0 M ) and A(w, MjOLS ) were calculated. For M (in A(w,9m) 
and A(w,9m,ols)), we used a discrete uniform distribution on the grid 
Wj = (wn, 1 — Wi 2 ) {wn = 0.05 ■ i, i = 1,2, ...,19). As expected, the 
naive nonparametric estimator A turned out to be clearly inferior to all 
other methods. For instance, for 9 = (0.1,0.1) and n = 50, the integrated 
mean squared error IMSE nonp = f A ^ E[(A(w) — A(w)) 2 ]dw is 74 times 
larger than IMSE nonp>OLS = E[(A OLS (w) - A(w)) 2 ]rfw, and IMSE p&r = 
f A E[(A(w, 9 M ) — A(w)) 2 ]dw is almost 7 times larger than the correspond- 
ing quantity (denoted by IMSE paurt oLs) for A(w, 9m,ols)- For larger sample 
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sizes the ratios 

' nonp 

IMSE nonp / IMSE noQPt OLS 

and 

'"par 

IMS E pax I IMS -Epar,OL5 

stabilize around the values of 64 and 30 respectively. Moreover, even if 
we compare the nonparametric OLS, Aqlsi with the parametric estimator 
A(w, 9 M ), we obtain a ratio of I M S E pai / 1 M S E nonp0 LS ~ 26 for large sam- 
ple sizes. We may thus conclude that using a good initial nonparametric 
estimator for the parametric method is essential. Detailed results on A(w) 
and A(w,6m) are therefore omitted, and we focus solely on the compari- 
son between Aqls an d A(w, 9m,ols)- Figures^), b) and c) show the ratio 
r = IMSE paT)0 Ls/IMSE nonPi oLS f° r the three choices of 9 as a function of 
n. In all three cases, r stabilizes around a value below 1. The numerical 
values are given in table 1. As a function of w, the relative precision of A(w) 
compared to A(w,9m,ols) depends on w and the shape of A. This can be 
seen in figures [5k), b) and c), where simulated values of 



r l w 



E 


(A(w, 9 m ,ols) ~ A(w) 


2" 


E 


(Aols(w) - A(w)) 2 





are ploted as a function of wi, for different values of n. Figures |6j [7]and|8j 
with estimates of A for 50 series of length n = 25 (Fig. a,b) and n = 200 (Fig. 
c,d) respectively, illustrate a further problem with the nonparametric OLS. 
For small sample sizes, Aqls is often not exactly convex, which means that 
it is, with relatively high probability, not a proper dependence function. By 
definition, this problem does not occur for A(w,9m,ols)- Finally, boxplots 
of 9\ and §2 for the case with 9 = (0.05, 0.9) are given in figures 9^,) and 
b) respectively. One can see in particular that for small sample sizes the 
distributions of 9\ and #2 are skewed to the right and left respectively. This 
is due to 9\ and 9 2 being close to the border of the parameter space, and the 
restrictions 6> 1; 9 2 > and 9\ + 9% < 1. 

6 Final remarks 

In this paper we considered estimation of extreme value copulas based on 
parametric models that are defined in terms of the spectral measure. This 
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approach is very flexible, and in principle any type of dependence between 
extremes can be captured. The method is not restricted to the case where the 
marginal distributions are known, since any nonparametric estimator A can 
be used in the projection. Theorem 1 and Corollary 1 apply (with a (v, w) 
replaced by the corresponding asymptotic covariance function) whenever a 
functional limit theorem of the form given in ([8j holds for A. 

An important issue that would need to be addressed in future research is 
the extension to a larger class of copulas. In this paper, observations were 
assumed to be generated by an extreme value copula. In practice, an ex- 
treme value copula is usually reached only asymptotically (for multivariate 
maxima). In analogy to nonparametric extreme value copula estimators, con- 
sistent parametric methods will have to be developed for such situations. A 
further question is model choice, i.e. the question how to decide on the num- 
ber and type of spectral measures to be used as a basis. For data generated 
by an extreme value copula, standard methods such as AIC or BIC ([T].[l6]) 
may be useful. In the more general situation where an extreme value copula 
is only reached in the limit, the question is more complex. 
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n 


9 = (0.1,0.1) 


9 = (0.05,0.9) 


9 = (0.8,0.1) 


Zo 


U.o I 


U.5U4 


n /ion 


50 


0.444 


0.622 


0.652 


100 


0.512 


0.636 


0.696 


200 


0.613 


0.652 


0.767 


400 


0.756 


0.674 


0.815 


800 


0.840 


0.749 


0.876 


1600 


0.899 


0.805 


0.898 


3200 


0.903 


0.868 


0.900 


6400 


0.902 


0.887 


0.903 



Table 1: r = IMSE pai ,oLs/ 'IMSE mOL s for 9 = (0.1,0.1), (0.05,0.9) and 
9 = (0.8, 0.1) respectively, and sample sizes n = 25 • 2 j (j = 0, 1, 8). 



24 



(a) theta=(1,0) 



(b) theta=(0,1) 




Figure 2: Simulated samples X; = (Xn,X i2 ) ~ C(-,9) (i = 1,2,.. .,n) with 
n = 1000, and 9 = (1, 0), (0, 1), (0, 0) and (§, |) respectively. 
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(a) theta=(1 ,0) 



(b) theta=(0,1) 
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(a) IMSE-ratio as a function of n 
(theta=(0.1,0.1)) 
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(b) IMSE-ratio as a function of n 
(theta=(0.05,0.9)) 
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(c) IMSE-ratio as a function of n 
(theta=(0.8,0.1)) 




\ 1 1 1 r 

1000 2000 3000 4000 5000 6000 



Figure 4: Ratio of parametric and nonparametric IMSE for (a) 9 = (0.1, 0.1), 
(b) (0.05,0.9) and (c) (0.8,0.1) respectively. 
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(a) MSE-ratio as a function of w 
(theta=(0.1,0.1)) 



(b) IMSE-ratio as a function of n 
(theta=(0.05,0.9)) 





w 



Figure 5: Ratio of parametric and nonparametric MSE, r(w), for (a) 6 = 
(0.1,0.1), (b) (0.05,0.9) and (c) (0.8,0.1), and sample sizes n =25, 50, 100, 
200, 400, 800, 1600, 3200 and 6400 respectively. 
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(a) A and 50 nonparametric OLS-estimates 
(n=25, theta=(0.1,0.1)) 



(b) A and 50 parametric OLS-estimates 
(n=25, theta=(0.1,0.1)) 




;ure 6: 50 estimates A OL s(w) (fig. (a) and (c)) and A (w,9m,olsj f° r 
- (0.1, 0.1) and n G {25, 200}. The black line represents the true function 
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(a) A and 50 nonparametric OLS-estimates 
(n=25, theta=(0.05,0.9)) 



(b) A and 50 parametric OLS-estimates 
(n=25, theta=(0.05,0.9)) 




;ure 7: 50 estimates Aols(w) (fig- (a) and (c)) and A(w, 6m,ols) for 
= (0.05, 0.9) and n G {25, 200}. The black line represents the true function 
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(a) A and 50 nonparametric OLS-estimates 
(n=25, theta=(0.8,0.1)) 



(b) A and 50 parametric OLS-estimates 
(n=25, theta=(0.8,0.1)) 




0.0 0.4 0.8 0.0 0.4 0.8 



Figure 8: 50 estimates Aols'(w) (fig. (a) and (c)) and A(w, 9m,ols) for 
9 = (0.8, 0.1) and n e {25, 200}. The black line represents the true function 
A. 
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(a) Boxplots of estimated theta[1] lor n=25*(2"j) 
(theta0=0.05.0.9) 




(b) Boxplots of estimated theta[2] for n=25*(2"j) 
(theta0=0.05,0.9) 
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Figure 9: Boxplots of O ls,i (fig- a ) an d Qols,2 (fig- b) for the case with 
9 = (0.05, 0.9), and n = 25 ■ 2 J (j = 0, 1, 8). The horizontal line represents 
the true value of ^ (i — 1,2). 
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